Revivals, classical periodicity, and zitterbewegung of electron currents in monolayer 

graphene 
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Revivals of electric current in graphene in the presence of an external magnetic field are described. 
It is shown that when the electrons are prepared in the form of wave packets assuming a Gaussian 
population of only positive (or negative) energy Landau levels, the presence of the magnetic field 
induce revivals of the electron currents, besides the classical cyclotron motion. When the population 
comprises both positive and negative energy Landau levels, revivals of the electric current manifest 
simultaneously with zitterbewegung and the classical cyclotron motion. We relate the temporal 
scales of these three effects and discuss to what extent these results hold for real graphene samples. 
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I. INTRODUCTION 

Graphene, a single-atom layer of carbon atoms ar- 
ranged in a honeycomb lattice, has attracted enormous 
attention recently. It exhibits remarkable mechanic and 
thermal properties*, but the most studied aspect are 
its startling electronic properties. Graphene can carry 
high current densities, its resistivity is less than that of 
silver— (the lowest resistivy material at room tempera- 
ture), and when exposed to a magnetic field at suffi- 
ciently low temperatures it shows an anomalous quantum 
Hall effect 3 -. Many of these properties are attributed to 
graphene's peculiar band structure, quasi-free electrons 
in single-layer graphene propagating as masslcss Dirac 
particles^, which opens the possibility of probing rela- 
tivistic quantum mechanics effects in condensed matter 
systems. For example, perfect tunneling of electron wave 
packets through potential barriers, the Klein paradox, 
unobservable in particle physics, has been experimentally 
confirmed in graphene very recently^. Additionally, in- 
terference between positive and negative frequency com- 
ponents of wave packets superposition of both positive 
and negative energy eigenstates induces a highly oscil- 
latory motion, known as zitterbewegung (ZB), that ap- 
pears as damped, rapid oscillations of free Dirac parti- 
cles around their otherwise rectilinear average trajecto- 
ries. ZB has been the subject of intense analytical investi- 
gations, although direct observations remain elusive due 
to the enormous characteristic frequency and the small 
amplitude of this motion. However, ZB may influence 
quantities other than the average position of free Dirac 
particles^, e.g. the nontrivial behavior of the conductiv- 
ity in graphene^, and it has been argued that graphene 
in a magnetic field is a promising system for the experi- 
mental observation of ZB 11 . 

In this paper we describe how, depending on the elec- 
tron wave packet's initial conditions, electric current re- 



generation occurs in graphene under an external mag- 
netic field, due to revivals and fractional revivals of the 
wave function. In particular, it is shown that when the 
electrons are prepared in the form of wave packets as- 
suming a Gaussian population of only positive (or nega- 
tive) energy Landau levels, the presence of the magnetic 
field induce revivals of the electron currents, besides the 
classical cyclotron motion. When the population com- 
prises both positive and negative energy Landau levels, 
revivals of the electric current manifest simultaneously 
with ZB and the classical cyclotron motion. Electric cur- 
rents due to ZB of electrons in graphene have already 
been describe d 11 : 12 , but the occurrence of revivals, a phe- 
nomenon purely quantum mechanical in origin, has not 
been reported yet despite that these oscillations are slow 
enough and large enough to be detected. Next, we briefly 
review the phenomenon of wave packet revivals and then 
study how they show up in graphene. 

The time evolution of wave packets, relativistic and 
nonrelativistic, can be quite complex due to quantum 
interference. However, several types of periodicity may 
emerge depending on the character of the energy eigen- 
value spectrum. If the initial wave packet is a superpo- 
sition of eigenstates <p n (x) sharply peaked around some 
large central no, different time scales can then be identi- 
fied from the coefficients of the Taylor expansion of the 
energy spectrum E n around the energy E no 13 i 14: , 
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For instance, propagating wave packets initially evolve 
quasiclassically and oscillate with period Tqi = 
2ttH/\E' I, which is in accordance with the correspon- 
dence principle. Wave packets then spread and collapse, 
and the classical oscillations eventually damp out. At 
later times, multiples of the revival time Tr = 4ttH/\E£ |, 
collapsed wave packets (almost) regain their initial wave- 
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form and oscillate again with period Tq\. Moreover, at 
times that are rational fractions of Tr, wave packets split 
into a collection of scaled and reshifted copies called frac- 
tional revivals. Longer time scales can be defined be- 
yond Tr, v.g. the so-called superrevival time, at which a 
new cycle of full and fractional revivals commences again. 
Lastly, the term E„ in the expansion Eq. |T]) above just 
generates an unobservable overall phase (see, however, 
below). Revival phenomena have attracted considerable 
interest over the past decades. In particular, revivals and 
fractional revivals have been investigated theoretically in 
nonlinear quantum systems, atoms and molecules^, and 
observed experimentally in, among others, Rydberg wave 
packets in atoms and molecules, molecular vibrational 
states, and Bose-Einstein condensates 16 . Interestingly, 
methods for isotope separation^, wave packet control 18 , 
as well as for number factorization^ have been put for- 
ward that are based on revival phenomena. 
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FIG. 1: Time dependence of |A(t)| 2 for initial Gaussian wave 
packets with B = 10 T and (top panel) n = 15, a = 3, Tr ~ 
17 ps and (bottom panel) n = 11, a = 40, Tr ~ 11 ps. 



II. REVIVALS OF THE ELECTRIC CURRENT 
IN GRAPHENE 

To study wave packet revivals in graphene, we consider 
the Hamiltonian for electrons in graphene under a mag- 
netic field perpendicular to the plane and of intensity B. 
Following standard methods 11 , within the tight binding 
approximation and near the vicinity of the Ki point, one 
of the two inequivalent corners of the Brillouin zone, the 
Hamiltonian reads (we use units such that h = 1) 



Hi = v F 
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(2) 



where tt = k — eA is the quasiparticle momentum, A 
is the vector potential and e the electron charge. The 
electronic dispersion relation is given by E(k) w ±u^|k|, 
where k = (k x , k y ) is the momentum measured relatively 
to the Ki point and vf — 10 6 m/s is the Fermi ve- 
locity (se o 20 ' 21 for details). Upon introducing the mag- 
netic radius L = yJl/eB and using the Landau gauge 
A = (—By, 0,0), it is not difficult to show that the en- 
ergy eigenfunctions are given by 
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with n = 0, 1, • • • , s = ±1 for the conduction and va- 
lence bands, respectively, and k x E K. Here, /„(£) = 
e~^ 2 H n (0/VLC n , £ = y/L- k x L, C n = \/2™n!0F, 
and -ff n (£) are the Hermite polynomials. The energy 
spectrum is, in turn, E„ tS = sE n = sfly/n, with Q = 
\/2v F /L. 

We shall construct the initial wave packets as the linear 
combination 



(k x )<p k ^ n , s (x,y). (4) 



Two different situations will be considered. First, wave 
packets with s = 1, centered around given ko x and no 
with coefficients c ni s=-i{k x ) — and c n . s= i(k x ) — c n (k x ) 
Gaussianly distributed as 



Cn(k x ) — 



(5) 



The classical period and the revival time yield straight- 
forwardly Tq\ = A'k^Jtiq/VL and Tr — 167rrig^ 2 /i7, re- 
spectively. In this case, ZB is not possible because the 
wave packets contain positive energy states only. Sec- 
ondly, we shall consider a superposition state of two wave 
packets with opposite s and centered around the same 
ko x and no as before, with coefficients c ntS= ^i(k x ) = 
(k x ) Gaussianly distributed as given by 
([5]). In this second case, simultaneously with the classi- 
cal periodicity and revivals, ZB shows up with a period 

1/2 

Tzb = K/E no = 7r/f2n ' that is obtained from the first 
term in the expansion Eq. {!]), which no longer acts as 
an unobservable overall phase. 

To account for the regeneration of the time-evolving 
wave packets we shall use the autocorrelation function 23 

A(t) = J dxdy^*(x,y,0)^(x,y,t). (6) 

An alternative approach based on information entropies 
has been recently proposed 2 ^,. After defining 

f°° * n v n s 1 ("-o> 2 o) 2 
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one readily finds 



A{t)=Y J Un 
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The occurrence of revivals corresponds to the return of 
|A(i)| 2 to its initial value of unity, and the fractional 
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FIG. 2: Time dependence of electric currents j x (blue, dotted 
line) and j y (red, solid line) in graphene for B — 10 T, no = 
15, a — 3, and s = 1. (a) First classical periods of motion 
with To — 279 fs. (b) Long-time dependence with Tr ~ 17 
ps. The classical periods and the main fractional revivals are 
indicated by vertical dotted lines. 
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FIG. 4: Time dependence of electric current j y in graphene 
for no = 15, a = 3, and s — ±1. Rest of parameters as in Fig. 
2. The vertical dotted lines stand for (a) ZB of electrons with 
Tzb — 4.7 fs, (b) classical motion, and (c) revival behavior. 



the currents 



revivals to the appearance of relative maxima in |A(t)| 2 . 
This is shown in figure [1] for B — 10T and two different 
initial wave packets. In the top panel, where no = 15, 
a = 3, and s = 1, revivals can be clearly seen at Tp — 17 
ps and Tr/2, and fractional revivals at Tr/4 and 3Tr/4. 
In the bottom panel, however, upon changing to no = 11 
and a = 40, no wave packet regeneration takes place 
because of the smaller no and large a. 

To investigate the behavior of electron currents, we 
first compute the x and y components of the current. The 
electron velocity operators are given by Vj = i[H, rj]/H = 
vpOj, (j = x, y), where a x and a y are the Pauli matrices. 
When only positive (s — 1) or negative (s — —1) Landau 
energy levels are populated, upon expanding in the base 
^fc x ,ra,a; one finds for the expected temporal evolution of 
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FIG. 3: Same as in Fig. 2 but with no = 11 and a = 40. 
Tci ~ 239 fs and T R ~ 11 ps. 



j x = -eV x (t) = eSV F ^ Un-l,nCOS [{E n - E n -l)t] , 

71=1 
CO 

jy = -ev y (t) = ev F ^ Un-i,n sin \{E n - E n -i)t] ■ (9) 

n=l 

The main features of the time evolution of j x and j v are 
extracted by using the expansion Eq. (p}. One finds 
E n — E n -i w E' na + E[[ o {n — no) from which the scales 
Tci and Tr arise, in agreement with our previous dis- 
cussion. It is shown in figure 2a how at early times, for 
no = 15 and (7 = 3,^ and j v follow fairly well the clas- 
sical cyclotron motion with a period Tci ~ 279 fs. In 
a few periods the wave packet enters the collapse phase 
and the quasiclassical oscillatory behavior of the currents 
vanishes, only to emerge later at half the revival time 
(see Fig. 2b). A revival of the electric currents can be 
observed at a time about Tr ~ 17 ps. By contrast, for 
n = 11 and a — 40, the quasiclassical behavior van- 
ishes right away (Fig. 3a), and the wave packet does not 
regenerate at long times (Fig. 3b) due to insufficient lo- 
calization. Notice that for revivals to be neatly observed, 
there is no need to reach very large values of no as long 
as a is sufficiently decreased. 

We now consider an initial superposition state of two 
wave packets with no = 15, a = 3, s = 1, and s = — 1. It 
is straightforward to obtain 

j x (t) = 

oo 

j y (t) = 2ev F ^ U n -i.n{sm [(E n + £7„_i)f] 

71=1 

+ sin[(£„-£ n _ 1 )t]}.(10) 

In this case, three different types of oscillatory motion 
arise. At very short times (in the femtosecond scale), the 
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FIG. 5: Time dependence of electric current j y in graphene for 
no = 11, cr = 40, and s = ±1. Rest of parameters as in Fig. 
2. The vertical dotted lines stand for (a) ZB of electrons with 
Tzb — 5.4 fs, (b) classical motion, and (c) revival behavior. 



electronic current is affected by ZB, as shown in Fig. 4a. 
Using again the expansion Eq. ([1]) leads to E n + E n+ i w 
2E no , from which T ZB = 2^j2E nQ = Tr/iln^ 2 ~ 4.7 
fs time-scale arises, which coincides with the expression 
previously derived on more general grounds. At medium 
times (Fig. 4b), ZB oscillations can be seen superimposed 
on the clearly visible, quasiclassical ones. Similarly, full 
and fractional revivals can be identified in the picosecond 
scale (Fig. 4c). Observec„ iS= _i(fc a ) = that, contrar- 
ily to what one would expect, the fractional revivals are 
clearly seen, i.e. they are not blurred due to the presence 
of two wave packets rather than a single one. Turning to 
no = 11 and a = 40, ZB persists but the quasiclassical 
periodicity fades away after the first three classical pe- 
riods (Fig. 5a), and the revivals are completely missed 
(Figs. 5b and 5c). 

We now turn to examine to what extent the above 
predictions hold for real graphene samples. It is well 
known that, in addition to the Ki point, there is a second 
inequivalent point in the Brillouin zone, K2, where the 
valence and conduction bands meet. The Hamiltonian 
for electrons at the K2 point is given b y 24 i 25 



H 2 = v F 





-7Tt + Z7T. 



-ir x — nr l 




(11) 



and hence H% = —Hj. The eigenvalues are the same 
than those of Hi but the eigenfunctions are different, 



fn(0 
Sfn-l(0 



(12) 



and the quantum velocity is given by v x = —Vpa x and 
v y = vpOy. Electrons around the Ki and K2 points are 
excited to levels associated with cp and £, respectively. 
Therefore, the total electric current is the sum of two 
contributions coming from two different wave packets: 
'F as given by Eq. (01 and $ obtained after replacing 



FIG. 6: Influence of the broadening of the Landau levels for 
no = 15 and B = 10 T. Rest of parameters as in Fig. 2. For 
r = 0.7 meV < r max the revival behavior is still observable 
(panel (a)). For T = r max most of the revivals are missed 
(panel (b)), but revival behavior can still be appreciated in a 
logarithmic scale (panel (c)). 



9k x ,n,s by (,k x ,n,s in Eq. (j4|), the coefficients c n , s (k. x ) be- 
ing common to both wave packets (see 2 ^) because the 
energy spectra of both Hamiltonians are identical. The 
contribution to the electric current from these two wave 
packets is equal^i. 

Another feature of real graphene that might have an 
impact on revival phenomena is the existence of a small 
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gap, of the order of 10 meV, caused by impuritie. 
interactions, etc . 27 ' 28 . However, the periods Tbi and Tr, 
will be the same than those calculated in the gapless 
model due to the fact that they are derivatives of the en- 
ergies, and only the Tzb period will change quantitatively 
with a new interband frequency 0, = (fl 2 + E 2 ^) 1 / 2 /h—, 
A second consequence of the imperfections in graphene is 
the broadening of the Landau levels, which may strongly 
influence the revivals of the electric current calculated 
using idealized, delta-like Landau levels. Assuming finite 
widths for all energy levels, characterized by broaden- 
ing parameters r„, the Landau levels E n are replaced 
by complex energies E n = E n + iTrr^. The temporal 
evolution of the currents can then be calculated as 

00 

jx = eSVF^Un-LnCOsKEn-En^e-^'^-^, 

71=1 
OO 

j y = ev F ^2 t/ n -i,n sin [(E n - E n _x)t] 

n=l 



-(r„+r„_! 



)f!3) 



We have estimated the maximum width r max up to which 
revivals of the electric current are clearly observable for 
no = 15 and B = 10 T. Using the approximation r n = T 
for n around no, we obtain r max < 3.7 meV. The in- 
fluence of the broadening of the Landau level on the re- 
vivals of the electron current is illustrated in Fig. [6] For 
r = 0.7 meV < r max the revivals are still observable 
(panel (a)), but for T — r max most of them damp out 
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due to the exponential factor in Eq. (fT3f (panel (b)). 
The early time surviving revivals with T = r max are vis- 
ible when the current is plotted in a logarithmic scale 
(panel (c)). Therefore, the occurrance of revivals of the 
electron current depends critically on the broadening of 
the Landau levels. 



III. DISCUSSION AND CONCLUSIONS 

In conclusion, we have studied the dynamics of electron 
currents in graphene subject to a magnetic field. Several 
types of periodicity must be distinguished if the wave 
packets representing the electrons arc sufficiently local- 
ized around some large enough central quantum number 
riQ. In this case, currents initially evolve quasiclassically 
and oscillate with a period Tqi, but at later times the 
wave packet eventually spreads, leading to the collapse 
of the classical oscillations. At times that are multiple 
of Tr, or rational fractions of Tr, the wave packet (al- 
most) regains its initial form, and the electron current 
its initial amplitude. For this to occur, the presence of 
a quantizing magnetic field is necessary, for if B = 
the spectrum is continuous which rules out the possibil- 
ity of revivals^. Associated with the revival of the wave 
packet the quasiclassical oscillatory motion of the cur- 
rents resumes. Additionally, when both positive and neg- 
ative Landau levels arc populated, permanent ZB oscilla- 
tions are observed, in agreement with previous results^. 
Like Tqi and Tr, Tzb = Tt/E no can be simply obtained 
from the expansion Eq. (TT|), thus providing a unifying 



description of three main time scales for the study of 
electric currents in graphene. Interestingly, although all 
the several periods depend on B as 1/ \fB, their ratios 
Xci/TzB = Tr/Tci — 4no only depend on no. Moreover, 
since Tzb < Tqi < Tr current revivals are more accessi- 
ble to experimental probing. In particular, for the values 
used in this work Tr/Tzb = 16riQ ~ O(10 3 ). Interest- 
ingly, the analogy of Hamiltonian @ with the Jaynes- 
Cummings model allows us to relate ZB and Rabi oscilla- 
tions, as other authors have pointed out£ ' 29 i 30 . All these 
results carry over without change when a more refined 
modeling of graphene is used that accounts for the K2 
point of the Brillouin zone and the energy gap caused 
by inevitable imperfections and interactions. However, 
the broadening of the Landau levels destroys most of the 
revivals in the case of B = 10 T, for broadening param- 
eters r > r max = 3.7 meV. Since current experimental 
estimates of the line width are of the order of 5 meV for 
B = 10 T, the observation of revivals of the electric cur- 
rent in graphene would require cither stronger magnetic 
fields, which lead to lower Tr and line width^i, or the use 
of high quality graphene samples. We hope the present 
work will stimulate further experimental studies in this 
direction. 
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